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Abstract 

We study the tradeoff between reliability, data rate, and delay for half-duplex MIMO multihop 
networks that utilize the automatic -retransmission-request (ARQ) protocol both in the asymptotic high 
signal-to-noise ratio (SNR) regime and in the finite SNR regime. We propose novel ARQ protocol 
designs that optimize these tradeoffs. In particular, we first derive the diversity-multiplexing-delay 
tradeoff (DMDT) in the high SNR regime, where the delay is caused only by retransmissions. This 
asymptotic DMDT shows that the performance of an N node network is limited by the weakest three- 
node sub-network, and the performance of a three-node sub-network is determined by its weakest link, 
and, hence, the optimal ARQ protocol needs to equalize the performance on each link by allocating 
ARQ window sizes optimally. This equalization is captured through a novel Variable Block-Length 
(VBL) ARQ protocol that we propose, which achieves the optimal DMDT. 

We then consider the DMDT in the finite SNR regime, where the delay is caused by both the 
ARQ retransmissions and queueing. We characterize the finite SNR DMDT of the fixed ARQ protocol, 
when an end-to-end delay constraint is imposed, by deriving the probability of message error using an 
approach that couples the information outage analysis with the queueing network analysis. The exponent 
of the probability of deadline violation demonstrates that the system performance is again limited by 
the weakest three-node sub-network. The queueing delay changes the consideration for optimal ARQ 
design: more retransmissions reduce decoding error by lowering the information outage probability, but 
may also increase message drop rate due to delay deadline violations. Hence, the optimal ARQ should 
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balance link performance while avoiding significant delay. We find the optimal fixed ARQ protocol by 
solving an optimization problem that minimizes the message error subject to a delay constraint. 

Submitted to IEEE Trans. Info. Theory, April, 2011. 
I. Introduction 

Multihop relays are widely used for coverage extension in wireless networks when the di- 
rect link between the source and destination is weak. The coverage of relay networks can be 
further enhanced by equipping the source, relays and destination with multiple antennas and 
using multiple-input-multiple-out (MIMO) techniques for beamforming. Indeed, MIMO can be 
used either for beamforming, which improves the reliability, or for spatial multiplexing, which 
increases the data rate Q]|. These dual uses of MIMO gives rise to a diversity-multiplexing 
tradeoff in point-to-point and multihop MIMO systems, as discussed in more detail below. 

Recovery of packets received in error in multihop networks is usually achieved by automatic 
retransmission (ARQ) protocols. With an ARQ protocol, on each hop, the receiver feeds back 
to the transmitter a one-bit indicator signifying whether the message can be decoded or not. In 
case of failure the transmitter retransmits the same message (or incremental information, e.g., 
using a Raptor code [EJ [I3J ) until successful packet reception. The ARQ protocol can be viewed 
as either a one-bit feedback scheme from the receiver to the transmitter, or as a time diversity 
scheme employed by the transmitter. The ARQ protocol improves system reliability at a cost 
of increased delay. In order to design an effective ARQ protocol for multihop relay networks 
with MIMO nodes, first the fundamental tradeoffs between reliability, data rate, and delay of 
such systems must be determined, and then the protocol performance can be compared to this 
theoretical performance limit. 

A fundamental tradeoff in designing point-to-point MIMO systems is the tradeoff between reli- 
ability and data rate, characterized by the diversity-multiplexing tradeoff (DMT). The asymptotic 
DMT was introduced in Q focusing on the asymptotically high SNR regime. The finite SNR 
DMT was presented in [5]. The DMT has also been used to characterize the performance of 
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classical three-node relay networks, with a direct link between the source and the destination, 
when the nodes have single- antenna (SISO) or multiple antennas for various relaying strategies 
0, 0, 0. The DMTs for the amplify-and-forward (AF) and decode-and-forward (DF) relaying 
strategies are discussed in J6]|. Several extensions of the amplify-and-forward strategy have been 
proposed recently, including the rotate-and-forward relaying and flip-and-forward relaying 
IfTOl strategies, which employ a sequence of forwarding matrices to create an artificial time- 
varying channel within a single slow fading transmission block in order to achieve a higher 
diversity gain. A dynamic decode-and-forward (DDF) protocol, in which the relay listens to the 
source transmission until it can decode the message and then transmits jointly with the source, is 
proposed in ifTTTl and its DMT performance is shown to dominate the fixed AF and DF schemes. 
The DDF protocol is shown to achieve the optimal DMT performance in MIMO multihop relay 
networks in [8J. In this paper, we restrict our attention to multihop networks using the DF 
relaying strategy, since it enables us design an optimal ARQ protocol for MIMO multihop relay 
networks, as we will show later. 

Here we consider the diversity-multiplexing-delay tradeoff (DMDT), which was introduced 
in [fT2l as an extension of the DMT to include the delay dimension. Here the notion of delay 
is the time from the arrival of a message at the transmitter until the message is successfully 
decoded at the receiver, also known as the "sojourn time" in queueing systems. Delays are 
incurred for two reasons: (1) ARQ retransmissions: messages are retransmitted over each hop 
until correctly decoded at the corresponding receiver, and (2) queueing delay: ARQ results in a 
queue of messages to be retransmitted at the transmitter. Most works on DMDT assume infinite 
SNR for the asymptotic analysis, and the queueing delay has been largely neglected. This is 
because in the high SNR regime, retransmission is a rare event lfT3l . With this asymptotic infinite 
SNR assumption, presents the DMDT for a point-to-point MIMO system with ARQ, HH 
studies the DMDT for cooperative relay networks with ARQ and single- antenna nodes, and |fT51 
proves the DMDT-optimality of ARQ-DDF for the multiple access relay and the cooperative 
vector multiple access channels with single antenna nodes. However, the asymptotically high 
SNR regime does not capture the operating conditions of typical wireless systems in practice, 
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where errors during transmission attempts are not rare events |fT3l . Hence, to fully characterize 
the DMDT performance in the finite SNR regime, we must bring the queuing delay into the 
problem formulation. For a point-to-point MIMO system with a delay constraint and no feedback 
link, the tradeoff between the error caused by outage due to insufficient code length, and the 
error caused by delay exceeding a given deadline, has been studied in [fT6l using large deviation 
analysis. 

One of the goals of our paper is to study the effects of dynamic ARQ on the DMDT in relay 
networks. Hence, we consider a line network in which a node's transmission is only received by 
adjacent nodes in the line. This is a reasonable approximation for environments where received 
power falls off sharply with distance (i.e., the path loss exponent is large). For this multi- 
hop channel model we show that the optimal ARQ protocol requires dynamic allocation of 
the ARQ transmission rounds based on the instantaneous channel state, and we obtain its exact 
DMDT characterization. The more general case where non-adjacent nodes receive a given node's 
transmission is significantly more complicated, and the optimal DMT is unknown for this case 
even with a single relay [14J. 

The contribution of this paper is two-fold: (1) we characterize the DMDT of multihop MIMO 
relay networks in both the asymptotically high SNR regime and in the finite SNR regime where, 
in the latter, queuing delay is incorporated into the analysis; (2) we design the optimal ARQ 
protocol in both regimes. Our work extends the DMDT analysis of a point-to-point MIMO system 
presented in |fT3l to MIMO multihop relay networks. In the first part of the paper, we derive the 
DMDT in the asymptotic high SNR regime, where the delay is caused by retransmissions only. 
For a certain multiplexing gain, the diversity gain is found by studying the information outage 
probability. An information outage occurs when the receiver fails to decode the message within 
the maximum number of retransmission rounds allowed. Based on this formulation, for some 
multihop relay networks a closed-form expression for the DMDT can be identified, whereas 
for general multihop networks, determining the DMDT can be cast as an optimization problem 
that can be solved numerically. The DMDT of a general multi-hop network can be studied by 
decomposing the network into three-node sub-networks. Each three-node sub-network consists 
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of any three neighboring nodes in the network and the corresponding links between them. The 
asymptotic DMDT result shows that the performance of the multihop MIMO network, i.e., its 
DMDT, is determined by the three-node sub-network with the minimum DMDT. The DMDT of 
the three-node sub-network is again determined by its weakest link. Hence, the optimal ARQ 
protocol should balance the link DMDT performances on each hop by allocating ARQ window 
sizes among the hops. From this insight, we present an adaptive variable block-length (VBL) 
ARQ protocol and prove its DMDT optimality. 

Next, we study the DMDT in the finite SNR regime, in which the delay is caused by both 
retransmissions and queueing. We introduce an end-to-end delay constraint such that a message 
is dropped once its delay exceeds this constraint. We characterize the finite SNR DMDT by 
studying the probability of message error, which is dominated by two causes: the information 
outage event and the missing deadline event, when the block length is sufficiently long [|4l. 
Our approach couples the information-theoretical outage probability Q with queueing network 
analysis. In contrast to the analysis under asymptotically high SNR, this does not yield closed- 
form DMDT expressions; however, it leads to a practically more relevant ARQ protocol design. 
The end-to-end delay that takes the queueing delay into consideration introduces one more 
factor into the DMDT tradeoff and the associated optimal ARQ protocol design. Specifically, 
allocating more transmission rounds to a link may improve its diversity gain and, hence, lower 
the information outage probability; however, it also increases the queueing delay and, hence, may 
also increase the overall error probability as more messages are dropped due to the violation 
of the deadline. Thus, an optimal ARQ protocol in the finite SNR regime should balance these 
conflicting goals: our results will show that this leads to equalizing the DMDT performance of 
the links. 

We formulate the optimal ARQ protocol design as an optimization problem that minimizes 
the probability of message error under a given delay constraint. The end-to-end delay constraint 
requires us to take into account the message burstiness and queueing delays, which are known 
to be the main obstacles in merging the information-theoretical physical layer results with the 
network layer analysis |[P71l . We bridge this gap by modeling the MIMO multihop relay network 
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as a queueing network. However, unlike in traditional queueing network theory, e.g., [[T8l , |[T9ll , 
the multihop network with half-duplex relay nodes is not a standard queue tandem, because node 
i along the multihop queue tandem must wait to complete reception of the previous message by 
the node i + 2 before it can transmit to node i + 1 in the tandem. Another difference between our 
analysis and traditional queueing theory is that we study the amount of time a message waits 
in the queue (similar to |fT3lO rather than just the number of messages awaiting transmission. 
This poses a challenge because the distribution associated with these random delays is hard 
to obtain ||20l , unlike the distribution of the number of messages for which a product form 
solution is available [18]. In [12111 delay is studied by using a closed queue model and diffusion 
approximation. We derive the exponent of the deadline missing probability in our half-duplex 
multihop MIMO network by adapting the large deviation argument used in |]22l . The expression 
of the exponent again demonstrates that the system performance (in terms of the exponent) of 
multi-hop network with half-duplex relays is determined by the three-node sub-network with the 
minimum exponent. 

The remainder of this paper is organized as follows. Section [II] introduces the system model 
and the ARQ protocol. Section [HI] presents the asymptotic DMDT analysis for various ARQ 
protocols while proving the DMDT optimality of the VBL ARQ. Section [IV] presents the finite 
SNR DMDT with queueing delays, including some illustrative examples. Finally, Section [V] 
concludes the paper and discusses some future directions. 



II. System Model and ARQ Protocols 



A. Channel Model 



ARQ 



ARQ 




Fig. 1. MIMO multihop relay network with ARQ. The number of antennas on the ith node is Mi, and H 6 £ M ;+i xA/ i j s 
the channel matrix from node i to node (i + 1). 
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Consider an iV-node multihop MIMO relay network. Node 1 is the source, node N is the 
destination, while nodes 2 through N — 1 serve as relays. Node i has Mi antennas for i = 
1, • • • ,N. The system model is illustrated in Fig. [Q We denote this MIMO relay network as 
(Mi, M 2 , ■ • • , M N ). At the source, the message is encoded by a space-time encoder and mapped 
into a sequence of L matrices, {Xij 6 q m i*t : / = x ( . . . where T is the block length, 
i.e., the number of channel uses of each block, and L is the maximum number of end-to-end 
total ARQ rounds that can be used to transmit each message from the source to the destination. 
The rate of the space-time code is R. 

We define one ARQ round as the transmitter sending a whole block code of the message 
to the receiver. We assume that the relays use the DF protocol: node i, 2 < i < N — 1, 
decodes the message, and reencodes it with a space-time encoder into a sequence of L matrices 
{Xij G C MlXT :/ = !,•••, L}. The channel between node i and node (i + 1) is given by: 



where Y iy i e C Mi+lxT , i = 1, • • • , N — 1, is the received signal at node (i + 1) in the 
Zth ARQ round. Channels are assumed to be frequency non-selective, block Rayleigh fading 
and independent of each other, i.e., the entries of the channel matrices H^i E C Mi+lxMi 
are independent and identically distributed (i.i.d.) complex Gaussian with zero mean and unit 
variance. The additive noise terms W it i are also i.i.d. circularly symmetric complex Gaussian 
with zero mean and unit variance. The forward communication links and ARQ feedback links 
only exist between neighboring nodes. 

Other assumptions we have made for the channel model are as follows: 
(i) We consider half-duplex relays, that is, the relays cannot transmit and receive at the same 



(ii) We assume a short-term power constraint at each node for each block code, given by 



transpose. A long-term power constraint would allow us to adapt the transmit power and 
achieve power control gain, as we briefly discuss later in the paper. In the following results 




1 < / < L, 



(1) 



time. 



K{tr(X] t Xi j)} < Mi, Wi,l. Here E{-} denotes expectation, and f denotes the Hermitian 
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we assume a short-term power constraint in order to focus on the diversity gain achieved 
by the ARQ protocol. 

(iii) We consider both the long-term static channel model, in which H it i = Hi for all I, i.e. 
the channel state remains constant during all the ARQ rounds for a hop and is independent 
from hop to hop; and the short-term static channel model, where H it i are i.i.d. but not 
identical for the same /. The long-term static channel assumption is the worst-case in terms 
of the achievable diversity with a maximum of L ARQ rounds [fT2l . because there is no 
time diversity gain. The long-term static channel model may be suitable for modeling indoor 
applications such as Wi-Fi, while the short-term static channel model suits applications with 
higher mobility, such as outdoor cellular systems. 

B. Multihop ARQ Protocols 

Consider a family of multihop ARQ protocols, in which the following standard ARQ protocol 
is used over each hop. The receiver in each hop tries to decode the message after or during 
one round, depending on whether the synchronization is per-block based or per-channel-use 
based. Once it is able to decode the message, a one bit acknowledgement (ACK) is fed back 
to the transmitter that triggers the transmission of the next message. After one ARQ round, 
if the receiver cannot decode the message, a negative acknowledgement (NACK) is fed back 
to the transmitter. Then the transmitter sends the next block of the code that carries additional 
information for the same message. The retransmission over the iih. hop continues for a maximum 
number of rounds, called the ARQ window size. Once the ARQ window size is reached 
without successful decoding of the message, the message is discarded, causing an information 
outage. Then the next message is transmitted. The sum of the ARQ window sizes is upper 
bounded by L > 0, where 

N-l 

Y,Li<L. (2) 

4=1 

We consider several ARQ protocols with different ways to allocate the available ARQ windows 
among different hops: 

(i) A fixed ARQ protocol, which allocates a fixed ARQ window size of Lj for the transmitter 
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of node i, i = 1, • • • , N — 1 such that J^iLi 1 Li = L. 
(ii) An adaptive ARQ protocol, in which the allocation of the ARQ window size per hop is not 
fixed but adapted to the channel state. The transmitter of a node can keep retransmitting 
as long as the total ARQ window size of L has not been reached. We further consider two 
types of adaptive ARQs based on different synchronization levels: 

(1) Fixed-Block-Length (FBL) ARQ protocol: The synchronization is per-block based. The 
transmission of a message over each hop spans an integer number of ARQ rounds. 

(2) Variable-Block-Length (VBL) ARQ protocol: The synchronization is per-channel-use 
based. The receiver can send an ACK as soon as it can decode the message, and the 
transmitter starts transmitting a new message without waiting until the beginning of the 
next channel block. VBL has a finer time resolution than FBL and is more efficient in 
using the available channel block, at a cost of higher synchronization complexity. 

We assume that the ARQ feedback links has zero-delay and no error. 

III. Asymptotic DMDT 

We characterize the tradeoff among the data rate (measured by the multiplexing gain r), the 
reliability (measured by the diversity gain d), and the delay by the asymptotic DMDT of a system 
with ARQ. Following the framework of flU and lfl2l . we assume that the rate of transmission 
depends on the operating SNR p, and consider a family of space time codes with block rate 
R(p) scaling with the logarithm of SNR as 

R(p)=r log p. (3) 

A. Diversity Gain 

In the high SNR regime, the diversity gain is defined as the SNR exponent of the message 
error probability [4]. It is shown in [4] that the message error probability P e (p) is dominated 
by the information outage probability P ou t(p) when the block-length is sufficiently large. In 
the following we make this assumption. The information outage event is the event that the 
accumulated mutual information at the receiver within the allowed ARQ window size does not 
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meet the data rate of the message and, therefore the receiver cannot decode the message. Hence, 
the diversity gain for a family of codes is defined as: 

d(r) ^ - lim l0 f*^. (4) 

p^oo log p 

The DMT of an M x x M 2 MIMO system is denoted by S MlM2 \r) and defined as the 
supremum of the diversity gain d{r) over all families of codes. DMT of a point-to-point MIMO 
system is characterized in [|4]| by the following theorem: 

Theorem 1. For a sufficiently long block-length, the DMT d^ Ml ' M2 \r) is given by the piece-wise 
linear function connecting the points (r, (Mi — r)(M 2 — r)), for r = 0, • • • , min(Mi, M 2 ). 

B. Asymptotic DMDT 

To characterize the asymptotic DMDT for a multihop network in the high SNR regime, we 
need the following quantity. Assume that the channel inputs at both the source and the relays are 
Gaussian with identity covariance matrices. Define M* = min{Mj, M i+ i}, for i — 1, • • • , N — 1. 
For the long-term static channel, let A^i, • • • , X^m? be the nonzero eigenvalues of HiH\, for 
i = 1, • • • , N - 1. Suppose Xij = p~ a ^, for j = 1, ■ ■ ■ , M*, % = 1, • • • , JV - 1. At high SNR, 
we can approximate the channel capacities C^iifj) = logdet (i + -j^H iH\ \ as Cj(JTj) = 
log p Si ( ai > where 

M* 

S 4 (a<) = X> - ay) 4 ", (5) 
i=i 

(x) + = max{x,0}, and the vector cti = [a iy \ - ■ ■ a iy M*}- This Si(ot,j) plays an important role 
in the asymptotic DMDT analysis. The closer the SNR exponents «j/s are to unity, the closer 
the channel matrix is to being singular. Similarly, we can define {a\ A in the short-term static 



channel model and the corresponding matrix aij 6 IR M ** xL as [cti\k,i — ci 



id J 

i,k' 

Proofs for the asymptotic DMDT analysis rely on the notion of decoding time, which is the 
time at which the accumulated information reaches R(p). In the case of the short-term static 

'Here the exponential equality = is defined as f(p) = p c , if linip-joo '"^^ = c. The exponential inequalities < and > are 
defined similarly. 
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channel, for the FBL ARQ and other block-based ARQ protocols, the decoding time for the ith 
node is given by 

U = inf jt G Z+ : Y^CiiH^) > rlogpj = inf |t e Z+ : ^StfaJ) > r| , (6) 

where Z + denotes the set of positive integers. For the VBL ARQ and other non-block-based 
ARQ protocols, the decoding time is given by 

U = inf It el : + {t - [t\ )#(a}* J+1 ) > r| , (7) 

where [x\ denotes the largest integer smaller than x. Similarly we can define the decoding time 
for the long-term static channel model. We can view the accumulated mutual information as a 
random walk with random increments Si (a') > and stopping boundary r. 

In the following, we first state our results for the three-node network (Mi, M 2 , M 3 ), and then 
extend them to the general iV-node network. 

1) Long-Term Static Channel: The DMDT of the fixed ARQ protocol in the case of the 
long-term static channel is given by the following theorem: 

Theorem 2. With the long-term static channel assumption, the DMDT of the fixed ARQ protocol 
for a three-node MIMO multihop network with window sizes L\ and L 2 , Li G Z + , L\ + L2 < L, 
is given by: 



4 Ml ' M2 ' M3) (r, L l7 L 2 \L, + L 2 <L)= gin (l) } 



(8) 



Proof: See Appendix lAl 

Consistent with our intuition, © shows that the performance of a three-node network is limited 
by the weakest link. This implies that if there were no constraint for the Lj's to be integers, the 
optimal choice should equalize the diversity-multiplexing tradeoff of all the links, i.e., 

d (M u M 2 ) (l_)= d (M2,M 3 ) ( 1_\ (9) 



M 

With the integer constraint we choose the integer Lj's such that the minimum of (f- MuMi+1 ^> (j^j 
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for i = 1, 2 is maximized. 

The DMDT of the FBL ARQ protocol is a piece-wise linear function characterized by the 
following theorem: 

Theorem 3. With the long-term static channel assumption, the DMDT of the FBL ARQ protocol 
for a three-node MIMO multihop network is given by 

,<™> (r , - {«""•"" © + © } • <io) 

Proof: See Appendix El 

The DMDT of the VBL ARQ protocol cannot always be expressed in closed-form, but can 
be written as the solution of an optimization problem, as stated in the following theorem. 

Theorem 4. With the long-term static channel assumption, the DMDT of the VBL ARQ protocol 
for a three-node MIMO multihop network is given by 

d m L M M ^ L)= inf h{{a .. }l (ID 

(£*i, ct2)eo 

where h({aij}) = Y^f=i St=i(2j — 1 + — M i+ i\)aij. The set O is defined as 



0±{{CLI,CL 2 ) G 



> > a i>M * > 0, % = 1, 2, ' / Q V 7- r <t , 

Si(ai) + 62(0:2) LJ 

an J this is the optimal DMDT for a three-node network in the long-term static channel. 

Proof: See Appendix O 

Note that the DMDT of the VBL ARQ protocol in the three-node network, under the long- 
term static channel assumption, is similar to the DMT of DDF without ARQ given in [8], with 
proper scaling of the multiplexing gain. The optimization problem in (fTTI) can be shown to be 
convex using techniques similar to Theorem 3 in 

We have closed-form solutions for some specific cases where the optimization problem in 
(fTTI) has a simple form and can be solved analytically. For example, for a (M l5 1, M 3 ) network 
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(fTTT) becomes 

d ( y B f M3 \r,L)= inf + M 3 a 2jl 



Oil, 1,02,1 



subject to - 'r± - ■ — — < — 

(1 - ai,i)+ + (1 - a 2) i)+ L 



a iA >0, 2 = 1,2. (13) 

The DMDT for this case (and two other special cases) is given by the following corollary: 

Corollary. With the long-term static channel assumption the DMDT of the VBL ARQ protocol 

(1) for a (Mi, 1, M3) MIMO multihop network is given by 

, min{Mi,M 3 }i^^, < r < LI 2; 
0, otherwise. 

(2) for a (1, M, 1) MIMO multihop network is given by 

, M^fr, < r < L/2; 

dfex ,1) (r,i) = < 7 " ~ (15) 

0, otherwise. 

(3) /or a (2, 2, 2) MIMO multihop network is given by 

si^WM <r<L/2; 



7(2,2,2)/ n 



2-r/L 

tgfc, L/2<r< 2L/3; 

fcfl, 2L/3<r<L; 

0, otherwise. 



2) Short-Term Static Channel: The DMDT of the fixed and the FBL ARQ under the short- 
term static channel assumption are similar to those under the long-term static channel assumption, 
with additional scaling factors for DMDTs of each hop due to the time diversity gain. 

Theorem 5. With the short-term static channel assumption, the DMDT of the FBL ARQ protocol 
for a three-node MIMO multihop network is given by 

•^V. - ^gu* it fi) } • 
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Proof: See Appendix iDl 

Theorem 6. With the short-term static channel assumption, the DMDT of the VBLARQ protocol 
for a three-node MIMO multihop network is given by 

d%r M (r,L)= inf (16) 

(tti,aa)6S 

where h ({ckJj}) = X/i=i Sj=i X^=i (2j — 1 + |Mj — M i+ i\) a 1 ^, and the set Q is defined as 

Q 4 {( ai , a 2 ) G R M ? xL x R M % xL : 

<!>••• >a| iM » >0,VM, t! + t 2 >L}. (17) 

77ie defined in 0, depend on the cti's. This is the optimal DMDT for a three-node MIMO 
multihop network in the short-term static channel. 

Proof: See Appendix El 

C. DMDT of an N-Node Network and Optimality of VBL 

Next, we extend our DMDT results to general iV-node MIMO multihop networks. Note that 
in our model, since each transmitted signal is received only by the next node in the network, 
the transmission over the ith hop does not interfere with other transmissions. We will show 
the DMDT of this more general network is a minimization of the DMDTs of all its three- 
node sub-networks, due to half-duplexing and multihop diversity. The multihop diversity [H 
captures the fact that we allow simultaneous transmissions of multiple node pairs in half-duplex 
relay networks. For example, while node i is transmitting to node (i + 1), node (i + 2) can also 
transmit to node (z+3). This effect allows us to split a message into pieces, which are transmitted 
simultaneously in the network to increase the multiplexing gain. Using this rate- splitting scheme, 
we can prove the DMDT optimality of the VBL ARQ protocol. Due to their fixed block length, 
we are only able to provide upper and lower bounds for DMDTs of fixed ARQ and FBL ARQ 
in an iV-node network. 

Theorem 7. With the long-term or short-term static channel assumption, the DMDT of the VBL 
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ARQ for an N-node MIMO multihop network is given by 

d%T' MN \r,L) = . min d^f^ M ^\r )L)) (18) 

and this is the optimal DMDT for an N-node network. 

Proof: See Appendix IB ■ 
Theorem [7J says that the DMDT of an iV-node system is determined by the smallest DMDT 
of its three-node sub-networks. The minimization in Theorem [7J is over all possible three-node 
sub-networks instead of pairs of nodes, due to the half-duplex constraint: each low-rate piece of 
message has to wait for the previous piece to go through two hops before it can be transmitted. 
Theorem |7J also says that the VBL ARQ is the optimal ARQ protocol in the general multi-hop 
network. 

Theorem 8. With the long-term or short-term static channel assumption, the DMDT of fixed 
ARQ for an N-node network is lower bounded and upper bounded, respectively, by 

4 Ml '-' Miv) (r^i,---^iV-i)>. min ( ^p r , Li , L l+l 



Li + Li+i < L max i 
(19) 



and 



1=1, ■■■ ,N—2 



where L max = rnax,^ 2 {Li + L i+1 }. 

Proof: See Appendix [Gl ■ 

Theorem 9. With the long-term or short-term static channel assumption, the DMDT of the FBL 
ARQ for an N-node network is lower bounded and upper bounded, respectively, by 

d%T' MN \r,L) > . min &^ M ^\r,L - N), (21) 

1=1, ■■■ ,iV — 2 



and 



S F M B r MN \r,L)<. min d^ L M ^ M ^\r, L). (22) 

1=1."' ,N—2 
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Proof: See Appendix El ■ 
An intuitive explanation for the DMDT optimality of the VBL ARQ is as follows. Recall that 
U is the number of channel blocks, including retransmissions, needed to decode the message 
over the ith hop. For a three-node network, we can illustrate the information outage region in 
the region of t\ x £ 2 values as in Fig. [2] The outage region of the VBL ARQ is smaller than 
those of the fixed and the FBL ARQ. Due to its per-block based synchronization, the outage 
region boundary of the FBL ARQ is a piecewise approximation to that of the VBL ARQ. In the 
high SNR regime, we formalize the above intuition in the following corollary to Theorem HI 

t 2 



l\ 
l\ 






' (t 1+ t 2 =L),VBL 
(t 1+ t 2 =L),FBL j\ 

l\ 



^^^^ Outage region for fixed ARQ (t,>L, ort 2 > L 2 ) 
^^^^ Outage region for adaptive ARQ ( t,+t 2 > L) 

Fig. 2. Outage regions of the fixed ARQ, and two adaptive ARQs: the FBL and the VBL ARQs, with L\ + L2 = L. 

Corollary. With the long-term or short-term static channel assumption, for an N-node MIMO 
multihop network, the DMDT of the FBL ARQ converges to that of the VBL ARQ when L — > oo. 
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Proof: Using d2U> and (|22i when L -> oo, 

^™_ a dS5^'^(r,L) > dKr^Cr^) > i=i 5_ 2 ff ,+ll%!) M[l - iV/L]) 

i=l,...,JV-2 VBL V ' 

(23) 
■ 

D. Power Control Gain with Long Term Power Constraint 

With the long-term power constraint and channel state information at the transmitter (CSIT), 
we can employ a power control strategy to further improve diversity. Let the SNR in the /th 
round be p(l) = p 9 ®, where p is the average SNR, and g(l) is the function defining the power 
control strategy. In the high SNR regime, similar to © we can approximate channel capacities 
as Cj(JTj) = \ogp s '^ ai \ where S^aij) = Yl!j=i (^(0 — a i.j) + ■ Hence, with power control, all 
the asymptotic DMDT results in the previous sections hold with Si(cti) replaced by S'^cti). 

E. Examples for Asymptotic DMDT 

DMDT of (4,1 ,3) multi-hop system, long-term static, L = 4 
3 

2.5 

2 

§ 1.5 
1 

0.5 

°0 1 2 3 4 

r 

Fig. 3. The DMDT for a (4,1,3) multihop network with L = 4. 




In this section we show some illustrative examples for the asymptotic DMDT. We first consider 
the long-term static channel model. For a three-node (4, 1, 3) multihop network with L = 4, Fig. 
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DMDT of (2,2,2) multihop relay system, long-term static, L = 4 




r 



Fig. 4. The DMDT for a (2,2,2) multihop network with L = 4. 

[3] shows the DMDT of the fixed ARQ with L\ = L 2 = L/2, of the per-hop-performance- 
equalizing L\ and L 2 satisfying ©, as well as the DMDTs of the FBL and the VBL ARQs. 
Note that the DMDT of the VBL ARQ in Fig. [3] is the optimal DMDT for the (4, 1, 3) network. 
We also consider a (2, 2, 2) network, whose DMDTs are shown in Fig. |4j 

Fig. [5] presents the three-dimensional DMDT surface of the VBL and the FBL ARQs, re- 
spectively, for the (4, 1,3) multihop network. Note that as L increases, the diversity gain at a 
given r increases for both the FBL and the VBL ARQ protocols. Also note that the DMDT 
surface of the FBL ARQ is piecewise and that of the VBL ARQ is smooth due to their different 
synchronization levels. Fig. [6] illustrates the cross sections of the surfaces in Fig. [5] at L = 2 and 
L = 10, which demonstrates the convergence of the DMDTs proved in Theorem IIII-C1 

Next we consider the short-term static channel model. The DMDT of the (4,1,3) multihop 
network using the FBL ARQ is shown in Fig. |7J Note that the asymptotic DMDT in the short- 
term static channel model is not necessarily a multiple L of the corresponding DMDT in the 
long-term static channel model, which differs from the point-to-point MIMO channel [fT2l . where 
the asymptotic DMDT in the short-term static channel model is a multiple L of the corresponding 
DMDT in the long-term static channel model. 
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DMDT of (4,1 ,3) FBL DMDT of (4,1 ,3) VBL 



3 3 




Fig. 5. The three-dimensional DMDT surface for a (4,1,3) network, with the FBL ARQ (left) and the VBL ARQ (right). 

DMDT of (4,1,3) System, L = 2 DMDT of (4,1,3) System, L = 10 




Fig. 6. The slices of the DMDT surface in Figure [5] at L — 2 (left) and at L = 10 (right). 

DMDT of (4,1,3) FBL, L = 4 




12 3 4 

r 



Fig. 7. The DMDT for a (4, 1, 3) multihop network in the long-term static channel versus that in the short-term static channel. 
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IV. Finite SNR DMDT With Delay Constraint 

In analyzing the finite SNR DMDT, we add a practical end-to-end delay constraint: each 
message has to reach the destination before the deadline; otherwise it is discarded. We char- 
acterize the finite SNR DMDT by studying the probability of message error. With the delay 
constraint brought into the picture, the probability of message error has two components: the 
information outage probability and the deadline missing probability, which we analyze using 
the finite SNR DMT introduced in [5] and the queueing network analysis, respectively. In the 
finite SNR regime, the multiplexing gain is defined as: 

r A ^ (24) 

i og2 (i + M rP y 

where M r . is the number of antennas at the receiver. In the following we only consider the 
long-term static channel model and the fixed ARQ protocol. We first introduce the queueing 
network model. 

A. Queueing Network Model 

The messages enter the network at the source node and exit from the destination node, forming 
an open queue. Messages arrive at the source node as a Poisson process with a mean message 
inter-arrival time of A blocks. As in the previous sections, the unit of time is one block of the 
channel consisting of T channel uses. The end-to-end delay constraint is k blocks. Each node 
can be viewed as a service station transmitting (possibly with several retransmissions) a message 
to the next node. The time node i spends to successfully transmit a message to node (i + 1) is 
called the service time of the ith node, which depends on the channel state and is upper bounded 
by the ARQ window size Lj. The allocated ARQ window sizes satisfy J^i^i — k- 

As an approximation, we assume that the random service time at node i for each message is 
i.i.d. with an exponential distribution of mean (the actual service time has value distributed 

in the interval of [0, Li]). Here fi(Li), which we derive later, is the actual average service time of 
the ARQ process when the ARQ window size is Lj. With these assumptions we can treat each 
node as an M/M/l queue. This approximation makes the problem tractable and characterizes 
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the qualitative behavior of MIMO multihop networks. The messages enter the buffer and are 
processed based on the first-come-first- served (FCFS) rule. We assume + fi(L i+ i) < A, 

i = 1, • • • , iV — 2, so the queues are stable, i.e., the waiting time at a node does not grow 
unbounded. 

Burke's theorem (see lfT8l ) says that in an M/M/l queue with Poisson arrivals, the messages 
leave the server as a Poisson process. Hence messages arrive at each relay (and the destination 
node) as a Poisson process with rate (1 — Pi)/\, where pi is the probability that a message is 
dropped. When the SNR is reasonably high, we can assume the message drop probability is 
small and hence 1 — p { & 1. 

B. Probability of Message Error 

Denote the total queueing delay experienced by the nth message transmitted from the source 
to the destination as W n , and the number of transmissions needed by node % to transmit the 
nth message to node {% + 1) as t l n , if the transmitter can use any number of rounds. For the 
nth message, if it is not discarded due to information outage, the total "service time" is S n = 
J^iLi 1 mm {^n5 ^i} and the random end-to-end delay is D n = W n + S n . Recall for the fixed ARQ, 
the message is dropped once the number of retransmissions exceeds the ARQ window size of 
any hop, or the end-to-end delay exceeds the deadline. Hence, the message error probability of 
the nth message can be written as 

P e = P {Ulfit > L^} + P {n£Lf{4 < n {D n > k}} . (25) 

The first term in (1231 ) is the message outage probability: 

P out ({L t }\p) ^ P {u^On >Li}}, (26) 

which is identical for any message n since channels are i.i.d. The second term in (1231 ) is related 
to the deadline missing probability, and can be rewritten as 

P {n^R < Li] n {D n > k}} = [l - p ou ,({L,}| P )]P {D n > k\nl^ {t n < z,}}. (27) 
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Define the stationary deadline missing probability: 

Pdeadline {{Li}\ P , k) = lim P {D n > k^f {** < U } } . (28) 



In the following we derive (|26l) and (|28l) . 

1) Information Outage Probability: Since channels in different hops are independent, (|26l 
becomes 

JV-1 N-l 



P out ({L^lp) = J2 P { t n> L i} = J2 P i L ^i( H i) < r l0 S2(l + M l+lP )} 
i=l i=l 
AT-l 



8=1 



which is a sum of the per-hop outage probabilities P outy i {L%\p)- Using results in [[51 for point 
to-point MIMO, we have 

7 (M i+ i - M< + 2Z - 1, ^[(1 + M m p) 6 < - (1 + M i+1 p)M 



Air 



■Pout.1 {Li\p) = sup I [ 
where the set B{ is given by 



^■■■C)^t=i r(M i+1 - Mi + 21 - 1) 

(30) 



(3D 

7(771, x) = J"* £ m-1 e~*(it is the incomplete gamma function, and T(m) = (m— 1)! for a positive 
integer m. For orthogonal space-time block coding (OSTBC), we can derive a closed-form 
Pout,i(Li\p) using techniques similar to |[51 : 

Pout,osTBc({^}|p) = ^p|r 8 L i log a + ^-11^11^ <rlog 2 (l + M mP )|, (32) 

where ||A|| denotes the Frobenius norm of a matrix A, and the spatial code rate r s is equal 
to the average number of independent constellation symbols transmitted per channel use. For 
example, r s = 1 for the Alamouti space-time code [|23l . When Hi is Rayleigh distributed, its 
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Frobenius norm has the Gamma(l, Mj • M i+t ) distribution. Hence, (1321 becomes: 



Pout,osTBc({^}|p) = £ (Mt . Mm _i)^ (yK 1 + M mP)^ - ll.^i • M mJ • (33) 

2) Deadline Missing Probability: For a three-node network with half-duplex relay, there is 
only one queue at the source that incurs the queueing delay. For given r and p, we can derive 
the stationary deadline missing probability using a martingale argument: 

Theorem 10. For a half-duplex three-node MIMO multihop network with Poisson arrival of rate 
A and ARQ rounds L\ and L 2 , the probability that the end-to-end delay exceeds the deadline k 
is given by 

P deadlme {{U}\ P , k) = MLl) + ML2) e-K^I)W-i). (34) 
Proof: See Appendix |U 

For general multihop networks with any number of half-duplex relays, the analysis for (|28l) is 
more involved. Due to half-duplexing the neighboring links cannot operate simultaneously, and 
this effect is not captured in the standard queueing network analysis. Here we adapt the proof 
in ll22ll . which uses large deviation techniques, to derive the following theorem for the exponent 
of the deadline missing probability in half-duplex relay networks: 

Theorem 11. For a half-duplex N-node MIMO multihop network, with Poisson arrival of rate A 
and ARQ rounds Li 's, the probability that the end-to-end delay exceeds the deadline k is given 
by 

lim lim \p{D n > kln^ 1 {£ < LM = -9*, (35) 

fc— >oo n— >oo k 

where 9* = mm 1 <i< N - 2 9i, and 9 { = + p,(L i+1 )] — 1/A, % = 1, • • • , N — 2. 

Proof: See Appendix |JJ 

This theorem again demonstrates that the performance of the iV-node multihop network with 
a half-duplex relay (here the performance metric is in terms of the deadline-missing probability 
exponent) is determined by the smallest exponent of each three-node sub-network. By Theorem 
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ITT1 for finite k we can approximate (|28l) as 



Pdeadline({L>i}\p, k) PS e 



(36) 



Also note that in the special case with iV = 3 nodes, for finite fc 



Pdeadline({Li}\p, k) ~ e 



M(£i)+M(I-2) A 



) 



(37) 



which is identical with ([341 up to a multiplicative constant [jti(Li) + /i(L 2 )]/A. This constant is 
typically not identifiable by large deviation techniques such as the one used in Theorem [TT] 

3) Mean Service Time Calculation: The above analysis requires p{Lj), which we will derive 
in this section. For a given t and message n, the cumulative distribution function (CDF) of t l n 
is given by 



where P out: i(t\p) is given in (1301) (or a term in the summation in (|33l) for OSTBC). Differentiating 
(1381) gives the desired probability density function (PDF): 



C. Optimal Fixed ARQ Design at Finite SNR 

Based on the above analysis, we formulate the optimal fixed ARQ design in the finite SNR 
regime as an optimization problem that allocates the total ARQ window size among hops to 
minimize the probability of message error subject to the queue stability and the end-to-end 



P K < t} = P {tdiHi) > rlog 2 (l + M i+lP )} = 1 - P mt ,i{t\p)i 



(38) 



where fi = — (1 + M i+ \p)~ — 1 . Using this we have 




(40) 
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delay constraint k: 



Him P out ({Li}\p) + [1 - P ut{{Li}\p)]Pdeadline{{Li}\p, k) 
{Li] 

N-1 (41) 

subject to 1 < fJ,(Li) < A, L, < k. 

i=l 



The terms in (1411) are given by (1291 ), (1331) . (1341) . and (|36l) . This optimization problem can be 
solved numerically. In particular, for a three-hop network with OSTBC, (I4TT) becomes 



{^} (MiM w - 1)! ' V P v 7 A 

N-1 

subject to 1 < p(Li) < A, L, < k. 

i=l 

As we demonstrate in the following examples, the information outage probability is decreasing 
in Li, and the deadline missing probability is increasing in Li. Hence the optimal ARQ window 
size allocation Lj at each node % should trade off these two terms. Moreover, the optimal ARQ 
window size allocation should equalize the performance of each link. 

D. Numerical Example 

We first consider a point-to-point (2, 2) MIMO system at the source and the destination. 
Assume p = 3 dB and A = 2 blocks. An OSTBC with r s = 1 is used, for which the information 
outage probability is given in (|33l) . The deadline constraint is k = 5 blocks. For r = 1, the 
information outage probability (|26|) and the deadline missing probability (1281) are shown in Fig. 
[8l Note that (|26|) decreases, while (|28l) increases with the ARQ window size. 

Next we consider the (4, 1, 3) MIMO multihop network. Assume p = 20 dB and A = 10 
blocks. An OSTBC with r s = 1 is used. The optimal fixed ARQ protocol is obtained by solving 
(|42|) . For r = 1 and a deadline constraint of k = 5 blocks, the optimal fixed ARQ has L\ — 2 
and L 2 = 3, and the optimal probability of message error is 0.1057. For r = 1 and k = 10 
blocks, the optimal fixed ARQ has L\ = 4 and L 2 = 6, and the optimal probability of message 
error is 0.0355. For all r and k, the probability of message error is plotted as a surface in Fig. 
|9l This surface is the DMDT for the three-node network in the finite SNR regime, which has 



26 



an interesting similarity to the high SNR asymptotic DMDT surface in Fig. [6l since indeed the 
high SNR DMDT represents the SNR exponent of the finite SNR DMDT. 




L 



Fig. 8. The information outage probability and the deadline missing probability for a (2, 2) point-to-point MIMO system, with 
r s = 1, r = 1, and SNR p = 3 dB, k = 5 blocks and A = 2 blocks. The minimum probability of message error is 0.4147 and 
is achieved by L — 2. 




Fig. 9. The probability of message error for the optimal fixed ARQ protocol in a (4, 1, 3) multihop network as a function of 
multiplexing gain r and delay constraint k; r s — 1, SNR p — 20 dB and A = 10 blocks. 



V. Conclusions 

We have analyzed the asymptotic diversity-multiplexing-delay tradeoff (DMDT) for the iV- 
node MIMO multihop relay network with ARQ, under both long-term and the short-term static 
channel assumptions. We show that the asymptotic DMDT can be cast into an optimization prob- 
lem that can be solved numerically in general, and closed-form asymptotic DMDT expressions 
are obtained in some special cases. We also proposed the VBL ARQ protocol which adapts 
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the ARQ window size among hops dynamically and proved that it achieves the optimal DMDT 
under both channel assumptions. We also show that the DMDT for general multihop networks 
with multiple half-duplex relays can be found by decomposing the network into three-node sub- 
networks such that each sub-network consists of three neighboring nodes and its corresponding 
two hops. The DMDT of the relay network is determined by the minimum of the DMDTs of the 
three-node sub-networks. We have also shown that the DMDT of the three-node subnetwork is 
determined by its weakest link. Hence, the optimal ARQ should equalize the link performances 
by properly allocating the per-hop ARQ window sizes dynamically. 

We then studied the DMDT in the finite SNR regime for fixed ARQ protocols. We introduced 
an end-to-end delay constraint such that a message is dropped once its delay exceeds the delay 
constraint. Since in the finite SNR regime retransmission is not a rare event, we incorporated 
the queueing delay into the system model, and modeled the system as a queueing network. 
The finite SNR DMDT is characterized by the probability of message error, which consists of 
the information outage probability and the deadline missing probability. While the information 
outage probability can be found through finite SNR DMDT analysis, we have also found the 
exponent for the deadline missing probability. Our result demonstrates that the performance 
of a multihop network with half-duplex relays in the finite SNR regime is also determined 
by the performance of the weakest three-node sub-network. It has been shown that, based on 
these analyses, the optimal fixed ARQ window size allocation can be solved numerically as an 
optimization problem, which should balance the per-hop diversity performance and avoid a long 
per-hop delay. 

The difficulty in merging the network layer analysis with the physical layer information 
theoretic results stems from the bursty nature of the source and the end-to-end delays. By 
modeling the multihop relay network with ARQ as a queuing network, we have tried to answer 
a question posed at the end of |fT3ll : how to couple the fundamental performance limits of 
general multihop networks with queueing delay. Our work provides a step towards bridging the 
gap between network theory and information theory. Future work includes developing an optimal 
dynamic ARQ protocol that can adapt to the channel state and the message arrival rate. The 
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problem can be formulated as a dynamic programming problem or analyzed using a heavy traffic 
approximation. 



Appendix A 
Proof of Theorem [2] 

With fixed-ARQ protocol and half-duplex relays, the system is in outage if any hop is in 
outage. The probability of message error P e (p), using the decoding time definition in ©, can 
be written as: 

P e (p) = P {{h >L 1 }U{t 2 > L 2 }} (43) 

= P{t 1 >L l } + P{t 2 >L 2 } (44) 

= P{L l S l {cc 1 ) < r} + P{L 2 S 2 {cx 2 ) < r} (45) 

p UJ (46) 



— mini 



(47) 



where (1441) is due to the independence of each link, and (|46l) follows from the method used in 
flU, and the fact that 

P {U >L l } = P |^ $(aj) < r J , (48) 

since «Si(a|) > and ^(a-) = 5j(o:j) for the long-term static channel. The last equality follows 
since when SNR is high, the dominating term is the one with the smaller SNR exponent. Using 
(|47l ) and the definition of diversity in © we obtain the DMDT stated in Theorem [21 

Appendix B 
Proof of Theorem [3] 

For the FBL ARQ protocol with two hops, the probability of message error is given by 

L-l 

P e {p) = P {h + t 2 > L} = J2? {h = k} P {t 2 > L - k} + P {tj > L} . (49) 

k=i 
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In the long-term static channel model we have 

P{h = k} = P{(k- 1)^(^0 < rlogp < fcdOEfi)} 

= p{c 1 (H 1 ) < ^ylogp} -P{Ci(i*i) < ^logp} (50) 

= p -^ (Mi ' M2) (^) - p - d(MM (i) = p ^ MiM2) (t^), 

which follows from the fact that S- Ml,M ' 2 \r) is monotone decreasing, i.e., S Ml ' M2 ^ (tzt) < 
d (M u M2) (ry If we plug (go) i nto (@9l) we get 

r log p 



P e (p) = P {d(fT > r log p} P <| C 2 (H 2 ) < y ! 



+ J2P{(L- k)C 2 {H 2 ) < r logp} Pit, = k) + P C^i) < I (51) 

fc=2 ^ ^ 

= p -min fc . 2 ,..., i _ 1 {d(M 1 ,M 2 )(^) id (M 1 ,M 2 ) fe)+ ,(M 2 ,M3)(^) id (M 2 ,M3)(_^)} 

= p -min il+i2 ^_ 1 ,^ {0 , 1 ,..., i _ 1} ,{d(M 1 ,M 2 )( ]2 _) +(i (M 2 ,M3,( t )^ 

where we have used the fact that S MuMi+1 \oo) = 0. From the definition of diversity in © the 
DMDT in Theorem [3] follows. 

Appendix C 
Proof of Theorem g] 

The decoding time of VBL ARQ is real, which differs from FBL ARQ. Since the long-term 
static channel has constant state, we can write the decoding time as t t = ■pjrf r \- Hence: 

(si i ) 

P e (p) = P{t 1 + t 2 >L} 

= P{(L- h)C 2 (H 2 ) < r logp < LC 1 {H 1 )} (53) 

r r 



= P{L< 



and dy^j^ 2 ' M3 \r, L) = inf{ a . d y e o h({ai t j}), where O is defined in (fT2]> . 

To prove that the DMDT of VBL ARQ is the optimal DMDT in an V-node network, we first 
provide an upper bound on the DMDT, and show that the DMDT of the VBL ARQ protocol 
achieves this upper bound. Our proof is for the short-term static channel in a three-node network, 



30 



as stated in Theorem [6j A similar (and simpler) proof can be written for the long-term static 
channel in a three-node network as stated in Theorem HI Assume that the source transmits for 
kT channel uses (kT < L) and the relay transmits in the remaining L — kT channel uses. Here 
k depends on the channel states and the multiplexing gain r. From the cut- set bound on the 
multihop network channel capacity, the instantaneous capacity of the MIMO ARQ channel is 
given by 



mm < max 

°x 1 ,,j=i,- ,lAI+i 



max 

^2 ,,!=!,■•■ ,L-W-1 



Y^I(X ltl ; Y lil \H 1 ,) + (k— [k\)I(X 1M1 ; Y 1M1 \H 1M1 ) 
1=1 

i-LfcJ-i 

J2 I(X 2X ,Y 2tl \H 2>l ) + (l-k+[k\)I(X 2 ^Y 2 ^ k \H 2>L _ k ) 



Since the capacity is maximized with Gaussian inputs, and linear scaling of the power constraint 
does not affect the high SNR analysis, the capacity C is upper bounded by 

{LfcJ L-[k\-l 
^2c l (H lil ) + (k-[k\)C 1 (H 1>k+l ), C 2 (H 2tl ) + (l-k+[k\)C 2 (H liL . 

1=1 i=i 

For any ARQ we can find a k* < L such that YlZi C 1 {H lil ) + (k*-[k*\)Ci{H 1)k * +1 ) = r log p. 
This means k* =t\. With this k*, the probability of message error is lower bounded by 

fz,-[*M-i 
rlogp> C 2 (H 2!l ) + (l-k*+[k*\)C 2 (H 1<L _ lk , l 



i=i 



l-L*iJ-i 



= P{<r> J2 5 2 (ay + (l-t 1 +L*iJ)S' 2 (a 



ng 



(55) 



i=i 



p ]{t 2 > L - ti} n g \ = P <{t! + 1 2 > L} n Q 



where Q = {(an, • • • , a N ^) G M. AI * xL x • • • x R A ^-i xL : a l iA >■■■> a{ M * > 0, Vz, /}. Hence, 
the diversity gain of any ARQ d( Ml ' M2 > Ms \r, L) of a three-node network is upper bounded by 



d^ M "' M3 \r,L)< inf ~h{a\,) 



(56) 
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with Q 2 = {h + t 2 > L} fl Q, which is the same as the set Q in (fTTT ), the DMDT expression for 
VBL ARQ in the short-term static channel. This shows that the DMDT upper bound is achieved 
by the VBL ARQ in the short-term static channel. This completes our proof. 

Appendix D 
Proof of Theorem [5] 

In the short-term static channel, for the FBL ARQ protocol with two hops (l50l) becomes 

k-1 "1 ( k 



P{ tl = k} = P ^d(/f M ) <r\ogp\-P lJ2 C ^ H ^i) < rlogp 
I i=i J I l=i 



Hence, the probability of message error can be written as 

P e (p) = P{C 1 {H ljl )>r\ogp}p\Y 4 C 2 {H %l )<T\og P 



1=2 



L-k 



+ J2 P \ S C 2 (H 2<l )<rlogp\P(t 1 = k) + Pl^2c 1 (H ltl )>rlogp 

k=2 U=fc+1 ) I 1=1 

By the definition of diversity, the DMDT in Theorem [5] follows. 

Appendix E 
Proof of Theorem [6] 

For a three-node network, we can break down the information outage event as a disjoint union 
of events that outage happens at the ith hop: 

P M ^^HA^' h{H ' } \ (57) 

where Qk — |X/t=i^ > Due to nonnegativity of U, Q\ C Q 2 . Hence, the minimization 
should be over Q 2 . Adding the ordering requirement on elements of {tt'j}, we have Theorem 
1 
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Appendix F 
Proof of Theorem [7] 

A. Upper bound 

We will first prove an upper bound for any ARQ protocol in an iV-hop network by considering 
a genie-aided scheme. For each % = 1, • • • , N — 2, consider the two consecutive hops from node 
i to node (i + 1) and then from node (i + 1) to node (i + 2). Assume a genie aided scheme 
where the messages are provided to node i, and the output of node (i + 2) is forwarded to the 
terminal node N. The maximum number of ARQ rounds that can be spent on this two-hop 
is L. The DMDT of this genie aided setup for any i, is an upper bound on the DMDT of the 
(Mi, • • • , M N ) system. The optimal DMDT of the {M h M i+1 , M i+2 ) system with L ARQ rounds 
is characterized in Theorem @J Hence, we have 

d^-^\r,L) < min i=li ... ^^'^(r,!), (58) 

where d^ Ml '- ' MN \r, L) is the DMDT of any ARQ protocol for an iV-hop network. 

B. The DMDT of VBL ARQ 

To be able to exploit the multi-hop diversity in the network, we use the following rate and 
ARQ round allocation scheme. First we split the original message of rate r log p into N/2 lower 
rate messages each having a rate of (r logp)/(iV/2) when iV is even (we split into (N — l)/2 
lower rate messages when N is odd). We pump these pieces of the original message into the 
network sequentially, and in equilibrium, they are transmitted simultaneously by adjacent pairs 
of nodes. 

Moreover, we require the number of blocks allowed for any two-hop transmission, from node 
i to node (i + 1) and then to node (i + 2), for all i = 1, ■ ■ ■ , iV - 2, to be L = L/(N/2) when N 
is even (or L = L/[(N — 1)/2] when iV is odd). This is equivalent to requiring the total number 
of blocks that each node i, i = 2, • • • , N, spends for listening and transmitting each piece of 
a message to be L. Note that with this constraint, the end-to-end total number of ARQ rounds 
used for transmitting each piece of the original message is upper bounded by L x N/2 = L 



33 

when N is even (or equals L x (N — l)/2 = L when N is odd). Hence, this scheme satisfies 
the constraint on the end-to-end total number of ARQ rounds. 

It is easy to see that the number of simultaneous transmission pairs we can have in an TV 
node network is N/2 when iV is even, and (N — l)/2 or (N + l)/2 when iV is odd. 

At the destination, all pieces of a message are combined to decode the original message. 
From the above analysis, the last piece of these low rate messages is received after at most 
L blocks, and the rate of the combined data is rlogp/(N/2) x (N/2) = rlogp when iV is 
even (similarly for odd N), which equals the original data rate. Hence this low rate message 
simultaneous transmission scheme meets both the data rate and end-to-end ARQ window size 
constraints. 

Now we study the outage probability P out (r) of this scheme. Define an outage event for any 
three-node sub-network consisting of nodes i, (i + 1), and (i + 2), for iV even, as: 

pi (r n A P f r/(iV/2) logp r /(N/2) log p L ] 
° uA ' ] I Q(H t ) + C l+1 (H l+1 ) > (N/2) J 



Ci(Hi) Ci+i(H i+ i) 



(59) 



and for N odd, similarly, as 



p (> L] ± p >/[(JV-l)/2]logp + r/[(iV-l)/2]logp ^ 



Ci(Hi) C i+1 (H i+1 ) \(N-l)/2] 

r log p r log p 



> L 



(60) 



Ci(Hi) Ci + \(H i+ i) 

Note that (|59l) and (|60l) say that by using this scheme, regardless of whether N is even or odd, 
the outage probability is as if we transmit the original message with data rate r log p over two 
hops with a total ARQ round constraint of L. From our earlier analysis for the VBL ARQ of a 
two-hop network, we have that as p — > oo, 

PUr,L)=p- d v BL (rj). (61) 



The system is in outage if there is an outage over any of the consecutive two-hop links from 
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the source to the destination. Using a union bound, we have 

N-2 

Pout(r,L)<J2PUr,L). (62) 

i=i 

As SNR goes to infinity, the right hand sum will be dominated by the slowest decaying term, 
which is the term with minimum d { ^f i+1 ' Mi+2 \r, L). Hence, the DMDT of this scheme is lower 
bounded by 

d Wu-,MN)^ L) > min d ^ L Mi+1 ' Mi+2 \r,L). (63) 

i=l,— ,N—2 

Together with the upper bound in (1581) . this shows that the presented scheme with the VBL 
ARQ achieves the optimal DMDT of an iV-hop network, and its DMDT is given by Theorem 

m 

Appendix G 
Proof of Theorem [8] 

The proof for the upper bound is identical to the one in Appendix IF-Al For the achievable 
DMDT of the fixed ARQ, we consider the following scheme with deterministic number of ARQ 
rounds: a node has to wait for at least Lj rounds over hop i for each piece of message, and 
we allow simultaneous transmissions to employ multihop diversity. Using this scheme, in steady 
state, the destination will receive one piece of the message every L max rounds (rather than L, if 
we do not employ multihop diversity). Now we divide the message into pieces with lower rates 
-^aa* r log p. Using this scheme, overall we will still achieve a rate of r log p in the steady state 
by transmitting these lower rate pieces. The outage probability of this scheme provides an upper 
bound on that of the fixed ARQ protocol: 

^Ll} \ Lmas f log q I 

P out (r,L ir -- ,L N ^) < J2 P i C{H ) >Ll }' (64) 

where \x] is the smallest integer larger than x. As SNR goes to infinity, the right hand sum will be 
dominated by the slowest decaying term, which is the term with minimum 



35 



and, hence, 



(Mi,-,Mn) 



V, L u ■ ■ • , L N _ X ) > min S M - M ^ ■ f ) 

i=l, ■■■ ,N—1 \ L Li J 



min d 

i=l ... N-2 



(Mi,Mi+l,Mi+2) I L mSLX 



: r T T 



Li + L i+ i < L n 



(65) 



which completes our proof. 



Appendix H 
Proof of Theorem [9] 

The proof for the upper bound is identical to Appendix IF-AI For the achievable DMDT of 
the fixed ARQ, again consider the same rate- splitting scheme in Appendix IF-Bl The difference 
here is that the number of ARQ rounds used is rounded up to be integer. For iV even, the outage 
probability can be written as 



PUr,L) = P 



r/(N/2) log(p) 
CAHA 



+ 



r/(iV/2)log(p) 



> L 



p , r/(N/2)log(p) + 1 + r/(JV/2)log(p) +1>L 

Ci(Hi) C i+ i(H i+ i) 
p j >log(p) + r log(p) +N > L 



CAH; 



Ci+l(Hi + i) 



(66) 



Note that L > N since we need at least N hops. The system is in outage if any three-node 
sub-network in outage. Using the union bound, again we have 

N-2 N-2 



i=i 



d ( ^t h+1 ' Mi+2 \r,L-N) 



i=l 



and 



d 



(M 1 ,-,Mjv), 
FBL 1 



r, L) > min d^^r 

' ^ - i=l,- JV-2 VBL 



(Mi,Mi+i,Mi+2) , 



r,L — N). 



(67) 



(68) 
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Appendix I 
Proof of Theorem [TO] 

Theorem U can be proved using Theorem 7.4.1 of 11241 . The theorem views the random queueing 
delay of the nth message as a reflected random walk. The deadline missing probability can be 
interpreted as a boundary hitting probability of the random walk, which can be obtained via a 
standard martingale argument which will not be repeated here. Note that the service time in a 
half-duplex two-hop network has a mean service time of fi(Li) +/x(L 2 ) blocks (approximate the 
service time as exponentially distributed) when conditioned on the event n^ 1 -^* < Lj}. The 
mean message inter- arrival time is A blocks. Using the mean service time, the mean inter-arrival 
time, and the delay deadline constraint k in Theorem 7.4.1 of [24j, we obtain the statement in 
Theorem \W} 

Appendix J 
Proof of Theorem [TTI 

The following proof is adapted from the proof in E21 . where a conventional queue tandem 
is considered. The conventional queue tandem is equivalent to a full-duplex multihop network, 
where the transmission (service) of node i for a message has to wait for the transmission of the 
previous message from node i to node (i + 1). However, in our problem, we have a half-duplex 
scenario, in particular, the transmission (service) of node i for a message has to wait for the 
transmission of the previous message over node (i + 1) to node (i + 2). The half-duplex scenario 
leads to a different and more complex queueing dynamic that we will study more precisely in 
the following. 

For node i, i = 1 ■ ■ ■ N — 1, N > 3, let the random variable S l n denote the service time of the 
nth message at node i, and A l n be the inter-arrival time of the nth message at node i. Due to 
the half-duplex constraint, there are iV — 2 queues at the source and node i, i = 2, ■ ■ ■ , N — 2. 
After the completion of transmission of the previous message, the message will be transmitted 
from node i to node (i + 1) and to node (i + 2), for % = 1, ■ • • , N — 2. Because of this queueing 
dynamic, the waiting time of the nth message at node i, W^, satisfies the following form of 
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Lindley's recursion (see 11221 ): 

W£ = (W^ 1 + 5i_ 1 + 52+ l 1 -iO+ i = l,-~,N-2, (69) 

where (s) + = max(x, 0). The total time a message spent for transmission from node i to node 
i + 2 is given by the waiting time plus its own transmission time 

JT n =Wi + S* n + S*\ i = l,---,N-2. (70) 

Note there are overlaps in these transmission times D l n 's we defined above, so their sum provides 
an upper bound on the end-to-end delay of the nth message: 

AT-2 

D n <Y,Dl (VI) 

i=i 

Next we will write -D^'s in (|7TT) explicitly using a non-recursive expression. The arrival process 
at node i is the departure process from node (i — 1), which satisfies the recursion 112211 : 

A\ = A*' 1 + Dl 1 - i = l,--,N-2, (72) 

where A l n is a Poisson process with mean interarrival time A. A well-known result from queueing 
theory [|22l states the following: if the arrival and service processes satisfy the stability condition, 
that is, the mean inter-arrival time A is greater than the mean service time n(Li) + n(L i+ i) at 
each of the queues i = 1, • • • , TV — 2, then Lindley's recursion (|69l has the solution: 

^ = 5gj(4,»-i + o£i-i- T i+i.»)' i = h'--N-2, (73) 

where the partial sums are defined as a\ p = Y7k=i $1 an< ^ T i, P = Y%=i Hence, from (f70|) and 
([73l . we have 

Dl i = mp{al n + o*£-T i jt+ltn ), i = l,---N-2. (74) 
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One the other hand, from (1721) we have 



, j ^ L + D P-Dn, i< P +i, 

T l,n = \ ( 75 ) 



0, otherwise. 



Plugging (1751) into (1741) , we have 



D^max^ + ^-rJ^-D^ + D}- 1 ), i = 2, • • ■ , N - 2. (76) 

JiS: n 

Moving D % ~ x to the left-hand- side, we obtain the recursion relation 

IK, ■ IK, ' f^(°L ■ u • /; :. ; )- i = 2,-~,N-2. (77) 

Now from dTJ) we have D]r l = max, : -, (rrj \ L + a)^, - ^T-d+i^)- Plu gg in S this int0 
(1771 ) we have 

^ + ^ = ^ <n (4,n + 4tn + + 4*-^ - ^ = 2 ' • • • ' N - 2 - (78) 

Repeating this operation inductively by expanding r 1-1 , , ln , we obtain 



JV-2 

z — 7 ji<-<jiv-i=n 



'JV-2 



(79) 



_i=l 

Note that in the non-recursive expression (1791) , we split the interval [l,n] by increasingly 
ordered integers ji, • • • , jnv-i) = ^> an d the summations of random variables over these different 
intervals are mutually independent. This decomposition enables us to adopt a similar large 
deviation argument as in E2l . to estimate the exponent 9* in the form of P(J2i=i 2 — ^) ~ 
exp(— 6*k) for large k. Following a similar argument as in G21 by finding a condition under 
which the log-moment generating function is bounded for each independent sum-of-random- 
variables when n — > oo, we can show that 



1 f N ~ 2 

lim -P > D* > k I = - min ^, (80) 

fc-s>oo k \~( J i=l,-,N-2 

where 6 { = sup{6> > : A r (-0) + A (5l+5l+ i ) (6') < 0}, i = 1, • • • , JV - 2, and the log- 
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moment-generating functions for the inter-arrival time A T (9) = log(l — 9\)~ 1 , for the service 
time A( 5l+ 5 l +i)(6 ) ) = log(l — 9[ft(Li) + fi(L i+ i)]) (if we approximate the sum service time 
to be exponentially distributed). We can further solve that 9{ = + — A -1 , 

i = 2, ■ • • , N - 2. Because of CD), P(D n > k) < P(Y!?=i D n > k )> and hence tne exponent 
for the deadline missing probability is bounded by —9* < — min i= i „. 5 at_ 2 9%- 

Now we prove the lower bound. Note that the end-to-end delay D„ is greater than the delay 
in any three-node sub-network: D n > D % n , i = 1, • • • , N — 2. Using a similar argument, we can 
show that the exponent for the probability that the delay in a three-node sub-network exceeds k 
is given by — 6>j. Hence, we have 

-9* = lim - log P(D n >k)> lim \ \ogP(D l n >k) = -9 { . (81) 

k— s>oo k fe— >oo k 

Inequality (I8TT) still holds if we take the maximum over all i on the right-hand-side: 

-9* > max -9i = - min 9 { . (82) 

i=l,...,JV-l i=l,-,N-l 

This completes our proof. 
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